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摘 要 : 为 探讨 全 球 气 候 变化 背景 下 多 元 线性 混合 模型 (IsoSource) 和 贝 叶 斯 混合 模型 (MixSIAR ) 解 
析 宁 夏河 东沙 地 柠 条 绵 鸡 儿 (Caragana korshinskii) 水 分 利用 策略 差异 及 适用 性 ,利用 和 毛 氧 稳定 同位 


素 技术 ,结合 直 


接 对 比 法 .IsoSource 模型 和 MixSIAR 模型 对 比分 析 不 同 坡度 (6"、10"、16" 和 24?) 样 


地 本 条 绵 鸡 儿 在 生长 季 不 同时 期 对 各 潜在 水 源 的 利用 率 ,并 评估 两 种 模型 的 植物 水 分 溯源 效果 。 
结果 表明 : 柠 条 绵 鸡 儿 对 不 同 土 层 深度 土壤 水 的 利用 存在 明显 的 季节 性 差异 ,生长 季 初 期 , 随 着 坡 
度 的 增 大 , 柠 条 绵 鸡 儿 对 中 层 土壤 水 的 利用 率 呈 现 出 先 增 大 后 减 小 的 趋势 ;生长 季 中 期 , 随 坡 度 增 
大 , 柠 条 绵 鸡 儿 主 要 水 分 来 源 由 深层 土壤 转移 至 浅 层 土壤 ;生长 季 后 期 , 柠 条 绵 鸡 儿 主 要 水 分 来 源 


随 着 坡度 的 增 大 由 浅 层 土壤 转移 至 深层 土壤 。 基 于 


直接 对 比 法 的 定性 判断 结果 ,TsoSource 模型 和 


MixSIAR 模型 计算 坡度 6°、10°* 和 16° 样 地 柠 条 锦 鸡 儿 主 要 水 分 来 源 利用 率 的 适用 性 均 较 高 ;而 Mix- 
SIAR 模型 计算 坡度 24" 样 地 柠 条 绵 鸡 儿 主 要 的 水 分 来 源 以 及 其 贡献 率 具有 更 高 的 可 靠 性 。Iso- 
Source 模 型 更 适合 解析 较 小 坡度 (6" 和 16?) 样 地 柠 条 绵 鸡 儿 的 水 分 利用 策略 ;而 MixSIAR 模型 解析 
较 大 坡度 (10? 和 24?) 样 地 柠 条 绵 鸡 儿 水 分 利用 策略 的 适用 性 更 好 。 研 究 结果 可 为 我 国 干旱 区 植 


物 水 分 来 源 鉴别 方法 的 选择 提供 科学 参考 。 


X 键 idl: IsoSource 模 型 ; MixSIAR 模型 ; 氢 氧 稳定 同位 素 ; 水 分 利用 策略 ; 宁夏 河东 沙 地 


文章 编号 : 


全 球 变 暧 加 速 了 地 球 系统 的 水 循环 ,使 全 球 和 
区 域 降水 格局 发 生变 化 "。 宁 夏河 东沙 地 地 处 毛 乌 
素 沙 地 西南 缘 ,在 空间 上 是 我 国 北方 农 牧 交 错 区 最 
有 具 过 渡 性 地 域 特征 的 沙 地 ,也 是 黄土 高 原 向 哪 尔 多 
斯 台地 过 渡 的 地 市 ,属于 气候 变化 的 敏感 区 ”。 研 
究 干旱 区 典型 固沙 植物 的 水 分 利用 策略 ,对 评 佑 降 
水 格局 变化 背景 下 该 地 区 的 生态 可 持续 发 展 具有 
重要 意义 ”。 稳 定 同位 素 技术 在 定量 分 析 植 物 水 分 
利用 策略 的 相关 研究 中 得 到 广泛 应 用 *“。 研 究 表 
明 , 除 极 少数 盐 生 . 旱 生 和 湿地 植物 外 ,陆地 植物 根 
系 吸收 的 水 分 在 到 达 叶 片 或 未 栓 化 的 幼 嫩 校 条 之 
前 一 般 不 发 生 同 位 素 分 馏 ”, 且 各 潜在 水 源 之 间 的 
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氨 氧 稳定 同位 素 组 成 存在 显著 差异 。 因 此 ,通过 比 
对 植物 木质 部 水 与 潜在 水 源 的 稳定 同位 素 组 成 ,可 
有 效 辨 别 和 区 分 植物 水 分 来 源 *。 

多 元 线性 混合 模型 (IsoSource) 和 贝 叶 斯 混合 模 
型 (MixSIAR ) 被 广泛 应 用 于 阐明 干旱 区 植物 的 水 分 
利用 策略 ,但 两 种 方法 的 量化 结果 可 能 存在 一 定 差 
祥 明 等 在 西南 喀斯特 地 区 利用 Mix- 
STAR #2 Fil soSource # Ht (ERK (Kalopanax sep- 
temlobus ) , # # (Toona sinensis ) Fil 14 # (Platycarya 
strobilacea ) 的 水 分 来 源 时 ,发 现 MixSIAR 模型 计算 
植物 对 各 水 源 的 利用 比例 适用 性 高 于 IsoSource 模 
型 ; 张 宇 等 利用 IsoSource 模 型 和 MixSIR .MixSIAR 


& [9-10] 
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模型 研究 了 红 瑞 木 (Swida alba) 的 水 分 来 源 , 表 明 
MixSIR 模型 能 够 更 精确 地 给 出 红 瑞 木 主要 的 吸水 
层 位 以 及 对 不 同 深度 土 层 土壤 水 的 吸收 利用 比例 ; 
Wang 等 后 利用 IsoSourece 模型 和 贝 叶 斯 混合 模型 
(SIAR MixSIR MixSIAR ) 研 究 了 黄土 高 原 半 干 旱地 
KEJP] (Vitex negundo ) MRTE (Artemisia gmeli- 
nii) 、 白 刺 花 (Sophora davidii) 的 水 分 来 源 ,发 现 
SIAR „MixSIAR 模型 比 MixSIR 模型 具有 更 好 的 植物 
水 分 溯源 效果 。 部 分 学 者 利用 IsoSource 模型 对 典 
型 固沙 植物 的 水 分 来 源 和 水 分 利用 特征 进行 了 研 
究 ", 但 该 模型 只 是 基于 质量 平衡 进行 近代 算法 ， 
存在 未 考虑 随机 测量 误差 .纳入 潜在 水 源 中 稳定 同 
位 素 组 成 的 不 确定 性 等 缺点 ,使 得 IsoSource 模 型 的 
量化 结果 可 靠 性 降低 ,因此 有 必要 对 不 同 模型 的 水 
源 区 分 性 能 进行 对 比 研究 ,分析 其 对 解析 植物 水 分 
利用 策略 的 适用 性 。 

本 研究 以 宁夏 河东 沙 地 典型 人 工 固沙 植物 柠 
条 绢 鸡 儿 为 研究 对 象 ,利用 氢 氧 稳定 同位 素 技术 、 
直接 对 比 法 .IsoSource 模型 和 MixSIAR 模型 分 析 生 
长 季 不 同时 期 不同 坡度 柠 条 绢 鸡 儿 的 水 分 利用 策 
略 差异 ,评估 两 种 模型 计算 植物 水 分 来 源 的 适用 
性 ,为 阐释 全 球 气 候 变 化 背景 下 干旱 沙 地 人 工 固 沙 
植物 水 分 利用 策略 提供 依据 。 
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注 :该 图 基于 国家 测绘 地 理 信息 局 标准 地 图 服务 网 站 下 载 的 


1 研究 区 概况 


试验 地 位 于 宁夏 回族 自治 区 灵 武 市 白 龙 滩 国 
RBA SREP EX (106° 22’ 03" E, 38°03’ 34N, WIR 
1150 m) (Al 1) ,该 地 区 地 处 毛乌素 沙 地 西南 缘 , 属 
温带 大 陆 性 气候 。 多 年 平均 气温 6.2~8.8 C ,多 年 平 
均 降 水 量 198.9 mm, 多 年 平均 潜在 蒸发 量 1928.4 
mm, 降水 主要 集中 在 7 一 9 月 ”%。 人 工 林 树种 主要 
A RY AR BBS JL 22 i AN (Euonymus maackii) 、 沙 束 
(Elaeagnus angustifolia) ~ WI] i (Robinia pseudoaca- 
cia) ği Ht 4% (Populus alba) | Vb $4 Œ (Calligonum 
mongolicum ) 和 和 细 枝 山竹 子 (Corethrodendron scopari- 
um) 等 。 地 带 性 土壤 以 淡 栗 钙 土 、 棕 钙 土 和 风沙 土 


2 数据 与 方法 


2.1 样品 采集 与 测定 

本 研究 在 柠 条 绵 鸡 儿 人 工 林 内 自然 分 布 的 4 个 
坡度 (6°、10° .16* 和 24°) 进 行 研 究 , 每 个 坡度 分 别 布 
设 3 块 20 mx20 m 的 样 地 , 柠 条 绵 鸡 儿 样 地 基本 情 
况 和 土壤 基本 理化 性 质 分 别 见 表 1 和 表 2。2018 年 
5 一 10 月 每 月 分 别 采 集 一 次 植物 木质 部 和 土壤 样 


21.17 
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图 号 为 GS(2019)1823 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 


1 人 研究 区 位 置 示意 图 
Fig.1 Schematic diagram of the location of the study area 
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表 1 不 同 坡度 柠 条 绵 鸡 儿 人 工 林 基本 情况 


Tab.1 Basic conditions of Caragana korshinskii plantations on different slopes 


YRC) 树龄 /a 平均 株 高 /m 平均 冠 幅 /m 平均 地 径 /cm 株 行列 距 /m 

6 3~10 1.87+0.02 2.43x2.19 1.70+0.08 3x4 

10 3~10 1.54+0.12 2.46x1.95 1.30+0.23 3x4 

16 3~10 1.58+0.11 2.45x2.20 1.30+0.22 3x4 

24 3~10 1.15+0.07 2.40x1.95 0.80+0. 17 3x4 

表 2 土壤 基本 理化 性 质 
Tab. 2 Basic physical and chemical properties of soil 
IEC) 容重 /gcm” 全 所 /mg ke! 全 磷 /mg .kg” 全 钾 /mg'kg 速效 磷 /mg kg” PACA mg kg” 

6 1.54 0.03+0.01 0.20+0.01 86.6543.19 47.98+40.97 100.94+27.68 
10 1.58 0.07+0.03 0.24+0.02 73.72+10.27 13.17+6.66 106.75+37.18 
16 1.54 0.10+0.02 0.19+0.02 62.36+4.22 15.05+8.59 102.82+24.89 
24 1.53 0.05+0.02 0.17+0.03 60.51+2.91 7.5340.37 98.94+8.36 


品 。 在 4 个 坡度 样 地 内 分 别 选 取 树龄 . 冠 幅 相 似 、 地 
径 相近 、 生 长 状况 良好 的 3 棵 植株 ,用 枝 剪 分 别 从 植 
株 4 个 方向 剪 取 中 部 直径 约 0.53 cm, 长 约 3~5 cm 的 
非 绿 色 栓 化 小 枝 ,迅速 剥 除 外 皮 与 韧 皮 部 , 装 入 10 
mL 的 玻璃 样品 瓶 中 ,并 用 Parafilm 封口 膜 密封 后 置 
于 随身 携带 的 保温 箱 (0~5 % ) 中 冷藏 ,每 个 坡度 样 
地 均 采 集 3 瓶 样品 作为 重复 。 同 时 ,在 所 选 柠 条 纺 
鸡 儿 灌 从 下 方 地 热 相 对 均一 的 地 点 ,用土 外 按照 0~ 
10 cm 10~20 cm 20~40 cm 40~60 cm, 60~80 cm, 
80~120 cm 、120~160 cm 和 160~200 cm 分 层 采 集 土 
壤 样 品 , 每 层 采集 2 瓶 土 样 ,立即 密封 在 已 标记 好 的 
样品 袋 中 , 放 入 保温 箱 中 。 所 有 植物 木质 部 和 土壤 
样品 带 回 实验 室 后 放 入 冰柜 (-20 % ) 中 冷冻 保存 直 
到 水 分 抽 提 。 在 无 林 空旷 地 内 放置 3 个 标准 雨量 简 
(直径 20 cm 的 圆 桶 上 放置 漏斗 ,下 端 与 雨量 简 相 
连 ,漏斗 口 放置 一 个 乒乓 球 ,防止 降水 样品 蒸发 ) 收 
集 降水 量 大 于 等 于 0.1 mm 的 降水 样品 ,每 次 降水 结 
束 , 立 即将 3 个 标准 雨量 简 中 的 降水 样品 混合 后 转 
和 10mL 玻 璃 样品 瓶 中 ,并 用 Parafilm 封口 膜 密封 后 
放 入 保温 箱 中 带 回 实验 室 , 置 于 冰箱 (4 % ) 中 冷藏 
保存 待 测 。2018 年 5 一 10 月 , 共 收 集 到 降水 样品 34 
个 。 利 用 Vantage Pro 2 小 型 气象 站 (Davis Inc, 
USA) ,每 隔 15 min 对 大 气温 度 .降水 量 等 气象 数据 
进行 记录 。 

利用 全 自动 真空 冷凝 抽 提 系统 (LI-2100 ,北京 
理 加 联合 科技 有 限 公 司 ) 抽 提 土 壤 水 和 植物 木质 部 
水 。 所 有 水 样 经 0.45 hm 醋酸 纤维 素 滤 膜 过 滤 后 ， 


利用 基于 稳定 同位 素 红外 光谱 技术 的 液态 水 同位 
素 分 析 仪 (T-LWIA-45$-EP,ABB-Los Gatos Research, 
CA) 测 定 降水 土壤 水 .植物 木质 部 水 样品 的 5D 和 
3"0 组 成 ,3D FSO 的 测量 误差 分 别 低 于 +1%o 和 
+0.3%o。 测 定 结果 以 相对 于 维也纳 标准 平均 海洋 水 
(Vienna standard mean ocean water, VSMOW ) 的 千 分 
差 值 (%o) 表 示 : 


Rogan kh 
3D( 或 8*0) = sample 


ye 2 1000% CL) 
Ry-smow 


式 中 ;Re 和 Ryswow 分 别 为 样品 和 维也纳 标准 平均 
海洋 水 中 的 稳定 同位 素 比 率 (D/H 和 ”*0/"*0)(%o)。 

利用 低温 真空 藻 馏 法 提取 植物 木质 部 水 时 ,可 
能 同时 藻 馏 出 与 水 分 子 具有 相似 光谱 吸收 峰 的 甲 
醇和 乙醇 等 有 机 物质 ,导致 稳定 同位 素 的 测量 误差 
超过 仪器 精度 """。 为 减 小 测量 误差 ,本 研究 利用 
Los Gatos 的 光谱 分 析 软 件 对 植物 木质 部 水 的 稳定 
同位 素 组 成 进行 校正 中 。 由 于 柠 条 绵 鸡 儿 是 干旱 
区 常见 的 旱 生 植物 ,而 旱 生 植物 的 根系 在 吸收 水 分 
的 过 程 中 会 发 生 氧 稳定 同位 素 分 馏 “ ,因此 ,本 研究 
采用 氧 稳定 同位 素 计 算 植 物 的 水 分 来 源 。 
2.2 研究 方法 

试验 地 地 下 水 埋 深 深度 为 18.0 m ,人工 固沙 植 
物 一 般 难 以 吸水 和 利用 地 下 水 ,因此 ,本 研究 只 将 
不 同 土 层 深 度 的 土壤 水 作为 潜在 水 源 。 

(1) 多 元 线性 混合 模型 (IsoSource ) 

将 柠 条 绵 鸡 儿 木质 部 水 8*0 值 作为 混合 物 ,不 
同 土 层 深度 的 土壤 水 及 其 8*0 值 作为 来 源 代入 模 
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型 ,将 来 源 增 量 设 为 1% ,将 质量 平衡 公差 设 为 
0.01 %o ,得 到 各 潜在 水 源 的 平均 水 分 贡献 率 。 

(2) 贝 叶 斯 混合 模型 (MixSIAR ) 

将 柠 条 绵 鸡 儿 木 质 部 水 8"*0 值 作 为 混合 物 数 
据 , 各 层 土 壤 水 8"0 值 作为 源 数据 ,并 将 分 馏 数 据 设 
为 0, 默 认 不 发 生 氧 稳定 同位 素 分 馏 , Error structure 
VEFE Resid Only ,马尔 卡 夫 链 蒙特 卡 罗 (MCMC ) 的 运 
行 长 度 设 为 normal, 模 型 结果 使 用 Gelman-Rubin 和 
Geweke 诊断 确定 MixSIAR 模型 是 否 收敛 ,如 果 Mix- 
SIAR 模型 不 收敛 则 增加 MCMC 的 运行 时 长 ,选择 
long 或 者 very long, 模 型 输出 结果 为 各 潜在 水 源 对 
应 的 平均 贡献 率 "。 

由 于 目前 植物 对 不 同 水 源 的 实际 利用 值 无 法 
直接 观测 Y, 因 此 ,本 研究 通过 评估 植物 木质 部 水 
氧 稳定 同位 素 的 观测 值 和 预测 值 之 间 的 匹配 程度 来 
间接 评估 IsoSource 模 型 和 MixSIAR 模型 的 性 能 3。 
将 测定 的 柠 条 强 鸡 儿 木 质 部 水 8"0 值 作 为 观测 值 
(O;) ,预测 值 LP) 的 计算 基于 柠 条 绢 鸡 儿 木质 部 水 
80 值 是 各 潜在 水 源 混合 物 的 假设 ”。 预 测 值 己 计 
算 如 下 : 


pi fd, (2) 


式 中 :n 为 植物 木质 部 水 的 潜在 水 源 个 数 (n=3);f 为 
IsoSource 模 型 和 MixSIAR 模型 计算 植物 对 第 i 个 洪 


(2) 模拟 优 度 指 数 (C): 


DL: z 0) 


C=1- 三 -一 一 (4) 
(0-0) 
i=1 
式 中 :0 为 观测 值 的 平均 值 。6G 值 越 接近 1 表示 模 
型 预测 的 可 靠 性 越 高 ;6 为 1 表示 模型 预测 非常 准 
确 ;G 值 为 负 值 表示 模型 预测 可 靠 性 不 高 。 
(3) 平均 绝对 百分比 误差 (MAPE): 
P,-0, 
0; 
MAPE<10% 表 示 高 度 拟 合 ;10%<MAPE<20% 表 
示 和 良好 拟 合 ;20%<MAPE<50% 表 示 可 行 拟 合 ; 
MAPE>50% 表 示 不 可 行 拟 合 。 
(4) 最 大 预测 误差 (MaxE): 
MaxE = max(P,- 0), i=1,2,.,n (6) 
(5) 最 小 预测 误差 (MinE ) : 
MinE = min(P,— 0), i=1, 2, .7 (7) 


3 结果 与 分 析 


3.1 降水 分 布 及 降水 中 5D 和 5*O 变化 特征 
试验 地 2018 年 5—10 月 植物 生长 季 总 降水 量 

为 172.1 mm, 主要 集中 在 7 一 9 月 (135.1 mm)( 图 2)， 

占 生长 季 总 降水 量 的 78.5%。 降 水 中 3D、5*0 的 变 


(5) 


-15 
MAPE = 1$ 


在 水 源 的 利用 率 (%);6, 为 各 潜在 水 源 的 氧 稳定 同 
位 素 值 (%o)。 

本 研究 根据 不 同 土 层 深度 土壤 水 5*0 值 的 相似 
性 及 变化 趋势 合并 相 邻 士 层 2 ,将 0~200 cm 的 土壤 
剖面 自 上 而 下 依次 划分 为 浅 层 (0~40 em) .中 层 
(40~120 cm) 和 深层 (120~200 em) 土壤 ,各 土 层 土壤 
水 8*0 值 取 对 应 深度 土壤 水 8*0 值 的 几何 平均 值 。 
本 人 研究 采用 Microsoft Excel 2019 软件 进行 数据 整 
理 ,在 SPSS 26.0 中 对 不 同 坡度 样 地 土壤 水 和 柠 条 强 
鸡 儿 木质 部 水 8*0 值 的 差异 显著 性 进行 单 因素 方差 
分 析 。 利 用 Origin 2018 软件 进行 图 像 绘制 ,图 表 中 
数据 为 平均 值 + 标准 差 。 
2.3 模型 评价 

(1) 均 方 根 误差 (RMSE ) : 


(3) 


式 中 :n 为 验证 样本 的 数量 ;0; 和 PP; 分 别 为 植物 木质 
部 水 530 的 观测 值 和 预测 值 。RMSE 越 低 ,表明 该 
模型 误差 越 小 ,预测 结果 越 准确 。 


化 范围 分 别 介 于 -74.1%o~7.4%o .-10.5%0~0.7%0Z [ia] 
(图 2)。 对 宁夏 河东 沙 地 事件 尺度 降水 中 5D、5"*O 
进行 线性 回归 分 析 , 拟 合 得 到 局 地 大 气 降 水 线 (LM- 
WL:8D=5.128*0-12.53, R=0.86,n=34)( 图 3)。LM- 
WL 的 斜率 和 截 距 均 小 于 全 球 大 气 降 水 线 (GMWL: 
8D=88”*O+10)”” ,说 明 雨 滴 在 非 饱 和 大 气 降落 过 
程 中 受到 了 较 强 烈 云 下 二 次 蒸发 效应 的 影响 ,降水 
中 稳定 同位 素 发 生 不 平衡 分 饱和 重 同 位 素 富 集 , 使 
48 LMWL WR mD 
3.2 不 同 坡度 柠 条 绵 鸡 儿 木质 部 水 6*O 变化 特征 
不 同 坡度 样 地 柠 条 绵 鸡 儿 木 质 部 水 8*0 差 异 显 
著 (P<0.05), 其 中 ,坡度 10°、16°* 和 坡度 24° 样 地 柠 条 
绵 鸡 儿 木 质 部 水 8*0 随 坡度 的 增 大 逐渐 减 小 (图 4)， 
变化 范围 分 别 介 于 -6.9%o~-2.1%o .-7.7%0~—-5.3%o 
和 -5.9%o~-0.49%o ,随时 间 的 变化 趋势 基本 一 致 , 生 
长 季 初 期 (5 月 ) 值 最 大 ,生长 季 后 期 (9 一 10 月 ) 值 最 
小 ;而 坡度 6° 样 地 柠 条 绵 鸡 儿 木质 部 水 80 的 变化 
范围 介 于 -7.8%o~-5.2%o ,和 后 长 季 后 期 值 最 大 ,生长 
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FHEWY 


35 7 (a) 气温 、 降 水 量变 化 


气温 /"C 


05-01 06-01 07-01 


一 气温 


降水 量 /mm 


09-01 10-01 11-01 


日 期 (月 -日 ) 


40 F (>) 降水 中 SD、8sO 变 化 


5D/%o 


06-01 06-21 07-11 07-31 


58°0/%o 


08-20 09-09 09-29 10-19 


日 期 (月 -日 ) 


图 2 不 同时 期 气温 、 降 水 量 及 降水 中 ;D .3"0 变化 


Fig.2 Variation of temperature, precipitation and ôD and 5"O in precipitation during different periods 


1207 e 降水 量 
100 上 ^ 坡度 6° F 
a 坡度 10。 一 LMWL: ôD=5.128!80-12.53 2 
80 | a HREIG --- GMWL: 8D=85"0+10 Po 
60 L * 坡度 24。 p A 
= 坡度 6 = 
40}+ «= 坡度 10° ik E 7 
x 20H © 坡度 16? A 3 
只 A aA 
0 上 
_20 L 
-40 L 
-60 L 
-80 L 
一 100 1 


1 1 1 
16 -14 -12 -10 -8 -6-4 20 2 4 6 8 10 
618O/%o 


注 :人 为 土壤 水 ; 口 为 柠 条 绵 鸡 儿 木 质 部 水 ;LMWL 为 局 地 大 气 降 
水 线 ;GMWL 为 全 球 大 气 降水 线 。 
图 3 局 地 大 气 降水 线 及 土壤 水 、 柠 条 锦 鸡 儿 木 质 
部 水 5D 和 56"*O 变化 


Fig. 3 Variation of local meteoric water 


line and ôD and 8"O in soil water and xylem 


water of Caragana korshinskii 


2 一 一 坡度 6。 一 一 坡度 16° 
So 一 e 一 坡度 10。 一 + 一 坡度 24° 
9 
% 2 
Z 
5o 
3 
6 
Š 
ip -8 

-10 L 1 1 1 1 
5 6 7 8 9 10 
月 份 


图 4 生长 季 不 同 坡度 森 条 绵 鸡 儿 木 质 部 水 6*O 变化 


Fig. 4 Variation of "O in xylem water of Caragana 


korshinskii at different slopes during the growing season 


季 中 期 (6 一 8 月 ) 值 最 小 。 
3.3 土壤 含水 量 与 土壤 水 6*O 变化 特征 

随 着 坡度 的 增 大 , 柠 条 绵 鸡 儿 样 地 土壤 含水 量 
逐渐 减 小 ,由 坡度 6° 的 (6.30+1.01)% 减 小 到 坡度 24? 
的 (1.33+0.58)%( 图 5)。 坡 度 24° 样 地 土壤 水 8*0 在 
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E HE: 宁夏 河东 沙 地 不 同 坡度 柠 条 绵 鸡 儿 (Caragana korshinskii) 水 分 利用 策略 差异 


土壤 含水 量 /% 土壤 含水 量 /% 土壤 含水 量 /% 
-2 0 2 4 6 8 1012141618 20 22 20 2 4 6 8 1012141618 20 22 202 46 81012141618 20 22 
0 0 Ore Ca ii 
Ai (a) 生长 季 初 期 O 生长 季 中 其 下 © 生长 季 后 其 
40 F 40 40 上 
60 F 60 60 H 
g 80 F § 80 § 80 F 
dy 100 上 i 100 a 100 上 
KK 120 + Kk 120 KE 120 上 
a 140 上 一 坡度 6。 5 140 一 坡度 6° 本 140 + 一 坡度 6° 
160 F 一 一 坡度 10 160 一 。 一 坡度 10 160 F —— 坡度 10。 
180 + 一 一 坡度 16 180 一 一 坡度 16 180 F 一 一 坡度 16。 
200 三 二 二 一 坡度 24 200 一 9 坡度 24 200 =y 坡度 24° 
220 + 220 220 = 


图 5 生长 季 不 同 坡度 土壤 含水 量变 化 


Fig. 5 Variation of soil water content at different slopes during the growing season 


0~40 cm Fil 40~120 em 土 层 最 高 ,120~200 cm 土 层 则 
表现 出 坡度 10° 最 高 .坡度 16° 最 低 。 不 同 坡度 柠 条 
绵 鸡 儿 样 地 土壤 水 8*0 随 士 层 深度 的 增加 逐渐 减 小 
(图 6)。 


27 一 s 坡度 6” 一 一 坡度 16° 
o L 一 e 一 坡度 10。 —— 坡度 24° 


5!80/%0 
上 
T 


月 份 
图 6 生长 季 不 同 坡度 土壤 水 5*O 变 化 


Fig. 6 Variation of soil water 6"O at different slopes 


during the growing season 


3.4 森 条 绵 鸡 儿 的 主要 吸水 层 位 

如 图 7 所 示 , 生 长 季 初 期 RE OOF FAT AR ANS 
儿 主 要 利用 中 层 (60~80 em) ,深层 (160~180 em) + 
IEK ,坡度 10?" 样 地 柠 条 绢 鸡 儿 主要 利用 中 层 (40~ 
60 cm 100~120 cm) 深层 (140~160 cm) 土 壤 水 ,而 
坡度 16° 样 地 柠 条 绵 鸡 儿 主 要 利用 中 层 (80~100 
cm、120 cm) 土 壤 水 ,坡度 24?" 样 地 柠 条 绢 鸡 儿 的 主 
要 吸水 层 位 在 60~80 cm 土 层 。 生 长 季 中 期 , 随 降水 
量 增加 ,坡度 6" 样 地 柠 条 绵 鸡 儿 与 土壤 水 存在 多 个 
交点 ,利用 直接 对 比 法 不 能 判断 出 水 分 来 源 ,坡度 
10? 样 地 柠 条 绢 鸡 儿 转 为 以 利用 浅 层 和 中 层 土壤 水 


和 深层 土壤 水 ,坡度 24° 样 地 柠 条 编 鸡 儿 增加 了 对 
浅 层 土壤 水 的 利用 ,同时 坡度 24° 样 地 柠 条 绵 鸡 儿 
木质 部 水 SO 与 中 层 土壤 水 存在 交点 ,表明 坡度 24? 
坡 样 地 柠 条 绵 鸡 儿 主 要 利用 浅 层 和 中 层 土壤 水 。 
生长 季 后 期 ,坡度 6" .10? 样 地 柠 条 绵 鸡 儿 转 为 以 利 
用 浅 层 (20~40 cm) 土 壤 水 为 主 , 而 坡度 16° 、24° 样 地 
柠 条 绵 鸡 儿 转 为 利用 更 深层 次 的 土壤 水 。 
3.5 IsoSource 模型 和 MixSIAR 模型 对 植物 水 分 来 
源 定量 分 析 

如 图 8 所 示 , 生 长 季 初 期 ,IsoSource 模 型 和 Mix- 
SIAR 模型 计算 坡度 6" 10 和 16? 样 地 柠 条 绵 鸡 儿 对 
土壤 水 利用 率 的 适用 性 均 较 高 , 随 坡 度 增 大 ,Iso- 
Source 模 型 计算 的 村 条 绵 鸡 儿 对 中 层 和 深层 土壤 水 
的 利用 率 分 别 高 于 MixSIAR 模型 的 计算 结果 
23.5% 、34.7% 和 24.4% , 而 MixSIAR 模型 计算 坡度 
24° 样 地 柠 条 纺 鸡 儿 土 壤 水 利用 率 的 适用 性 更 高 。 
生长 季 中 期 ,IsoSource 模型 和 MixSIAR 模型 计算 坡 
FE 10? 样 地 柠 条 绵 鸡 儿 对 浅 层 和 中 层 土壤 水 的 利用 
率 较 为 一 致 (为 91% ) ,但 IsoSource 模型 计算 坡度 
24? 样 地 柠 条 绵 鸡 儿 对 浅 层 和 中 层 土壤 水 的 总 利用 
率 高 于 MixSIAR 模型 的 计算 结果 7.6%。 生 长 季 后 
期 ,坡度 6? 样 地 柠 条 绢 鸡 儿 的 主要 水 源 由 深层 转移 
至 浅 层 土壤 ,IsoSource 模 型 计算 坡度 6°、16°? 、24° 样 
地 柠 条 编 鸡 儿 对 土壤 水 利用 率 的 适用 性 高 于 Mix- 
SIAR 模 型 ,坡度 16° 样 地 柠 条 绵 鸡 儿 对 深层 土壤 水 
的 利用 率 增加 至 最 大 (52.4% ) ,而 坡度 24° 样 地 柠 条 
绵 鸡 儿 的 主要 水 源 由 中 层 转 移 至 深层 土壤 
(73.9% ) ;坡度 10?" 样 地 柠 条 绢 鸡 儿 的 主要 水 源 由 中 


为 主 ,而 坡度 16" 样 地 柠 条 绢 鸡 儿 仍 主要 利用 中 层 


层 转 移 至 浅 层 十 壤 ,MixSIAR 模型 适用 性 更 高 , 柠 条 
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图 


7 不 同 坡度 本条 绵 鸡 儿 5"O 的 水 分 来 源 分 析 


x + 
titw 
一 一 土壤 水 ”一 一 木质 部 水 
as (a) 坡度 6° 生 长 季 初 期 220 - (b) 坡度 6" 生 长 季 中 期 220 - (c) 坡度 6° 生 长 季 后 期 
200 上 200 + 200 + 
180 + 180 + 180 + 
160 + 160 上 160 上 
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H 80+ H 80b # 80+} 
60 + 60 上 60 L 
40 + 40 + 40 上 
20 + 20 + t 20 b 
0 0 0 
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220  (d) 坡度 10° 生 长 季 初 其 220 「 (e) 坡度 10° 生 长 季 中 期 220 (f) 坡度 10° 生 长 季 后 其 
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Fig.7 Analysis of the water source of '°O in Caragana korshinskii at different slopes 
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高 ” 阳 等 ,宁夏 河东 沙 地 不 同 坡度 柠 条 绵 鸡 儿 (Caragana korshinskii) 水 分 利用 策略 差异 


= 0~40 cm mm 40~120 cm 
(b) IsoSource 模 型 -坡度 10。 


(a) IsoSource 模 型 -坡度 6 
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(c) IsoSource 模 型 -坡度 16?° (d) IsoSource 模 型 -坡度 24。 
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(g) MixSIAR 模 型 -坡度 16。 


初期 中 期 后 期 
生长 季 


(h) MixSIAR 模 型 -坡度 24° 
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图 8 IsoSource 模 型 和 MixSIAR 模型 计算 生长 季 不 同 坡度 柠 条 绢 鸡 儿 土壤 水 利用 率 


Fig. 8 IsoSource model and MixSIAR model for calculating soil water use in 


Caragana korshinskii at different slopes during the growing season 


纺 鸡 儿 对 浅 层 土壤 水 的 利用 率 为 38.3%。 
3.6 IsoSource 模型 和 MixSIAR 模型 的 植物 水 分 漳 
源 效果 

如 表 3 所 示 ,MixSIAR 模型 在 计算 坡度 10°、 
16° .24° 样 地 柠 条 编 鸡 儿 水 分 来 源 时 RMSE 更 小 ; 
IsoSource 模型 计算 坡度 6" 样 地 柠 条 绵 鸡 儿 水 分 来 
源 时 RMSE 更 小 。MixSIAR 模型 计算 的 坡度 6°、 
16° .24° 样 地 柠 条 纺 鸡 儿 的 6G 为 负 值 ,区 分 水 源 性 能 
可 靠 性 降低 ; 当 坡 度 为 24° 时 ,TsoSource 模型 的 
MAPE>50% 。IsoSource 模型 计算 坡度 6" 样 地 柠 条 


绵 鸡 儿 水 分 来 源 时 具有 更 好 的 植物 水 分 溯源 效果 ; 
MixSIAR 模型 计算 坡度 10° 样 地 柠 条 绵 鸡 儿 水 分 来 
源 的 济源 效果 更 好 ;计算 较 大 坡度 柠 条 绵 鸡 儿 水 分 
来 源 时 ,IsoSource 模 型 的 准确 性 低 于 MixSIAR 模型 。 


4 讨论 


坡度 \ 坡 向 等 地 形 条 件 对 地 表 光 、 热 .水 资源 具有 
青 分 配 作 用 ,进而 影响 林地 根系 的 生长 和 分 布 “”， 
而 植物 的 根系 分 布 特征 决定 植物 的 水 分 利用 策略 ”。 


表 3 IsoSource 模 型 和 MixSIAR 模型 的 预测 性 能 


Tab. 3 Predictive performance of IsoSource model and MixSIAR model 


评价 指标 、 fee 坡度 16 坡度 24 
IsoSource 模 型 MixSIAR 模 型 ”IsoSource 模 型 MixSIAR 模 型 ”J]soSource 模 型 MixSIAR 模 型 ”TsoSource 模 型 MixSIAR 模 型 
RMSE 0.00002 0.00157 0.00155 0.00087 0.00258 0.00002 0.00310 0.00004 
G 0.9997 -2.1632 0.1060 0.7645 0.9996 -8.2986 0.9998 =1.2755 
MAPE 0.0015 0.2036 0.0591 0.4094 0.0020 0.3396 1.3147 0.0245 
MaxE 0.0031 0.0007 0.0029 0.0016 0.0038 0.0017 0.0034 0.0013 
MinE —0.0017 —0.0007 —0.0017 —0.0016 0.0001 0.0001 —0.0062 0.0070 


HE: RMSE 为 均 方 根 误差 ;6 为 模拟 优 度 指 数 ;MAPE 为 平均 绝对 百分比 误差 ;MaxE 为 最 大 预测 误差 ;MinE 为 最 小 预测 误差 。 
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已 有 研究 “表明 , 随 着 坡度 的 增 大 , 细 根 根系 参数 
的 最 大 值 由 中 深层 土壤 向 浅 层 土壤 转移 ,呈现 出 表 
聚 现 象 ,并 且 随 士 层 次 度 的 增加 , 细 根 的 根 长 密度 、 
根 表 面积 密度 等 根系 参数 与 土壤 水 分 和 养分 的 变 
化 趋势 基本 一 致 , 呈 递减 趋势 。 在 整个 生长 季 , 标 
条 绢 鸡 儿 根据 土壤 含水 量 的 相对 高 低 ,在 不 同 土 层 
间 灵 活 转换 主要 水 分 来 源 。 生 长 季 初 期 PP ARS 
儿 对 中 层 土壤 水 利用 率 随 坡 度 增 大 呈 先 增 后 减 趋 
势 ; 生 长 季 中 期 , 森 条 绵 鸡 儿 的 主要 水 分 来 源 随 坡 
度 的 增 大 向 浅 层 土壤 转移 ,主要 利用 浅 层 和 中 层 土 
壤 水 ;生长 季 后 期 , 柠 条 绢 鸡 儿 主要 水 分 来 源 随 坡 
度 的 增 大 向 深层 土壤 转移 ,吸收 利用 深层 土壤 水 。 
这 可 能 与 柠 条 绵 鸡 儿 具 有 二 态 性 根系 结构 有 关 ”。 
朱 雅 娟 等 史 对 乌 兰 布 和 沙漠 村 条 锦 鸡 儿 的 研究 也 
表明 ,村 条 锅 鸡 儿 在 不 同 季 节 分 别 利用 浅 层 和 较 深 
层 的 土壤 水 分 ,推测 其 可 能 具有 双 型 根系 。 生 长 季 
初期 , 柠 条 绷 鸡 儿 需 水 量 少 , 水 分 来 源 单一 ; 而 生长 
季 中 期 ,一 方面 ,降水 补给 土壤 水 作用 明显 ,土壤 会 
水 量 显著 增加 ,植物 利用 浅 层 和 中 层 土壤 水 消耗 的 
能 量 相 对 较 小 ””, 男 一 方面 ,土壤 水 分 和 养分 含量 
供给 充足 , 细 根 数量 显著 增加 , 柠 条 绵 鸡 儿 优先 利 
用 浅 层 土壤 水 ” ;生长 季 后 期 , 柠 条 绢 鸡 儿 吸收 利 
用 深层 土壤 水 这 一 稳定 性 水 源 ,说 明 柠 条 绵 鸡 儿 采 
取 了 资源 依赖 型 的 水 分 利用 策略 ?”。 

IsoSource 模 型 和 MixSIAR 模型 在 量化 植物 水 分 
来 源 上 存在 一 定 差异 '"。IsoSource 模 型 和 MixSIAR 
模型 解析 坡度 6" .10" 和 16? 样 地 柠 条 纺 鸡 儿 主 要 水 
分 来 源 利用 比例 的 适用 性 较 高 。 然 而 这 两 种 模型 
在 计算 坡度 24° 样 地 柠 条 绵 鸡 儿 主 要 水 分 来 源 利 用 
比例 时 存在 一 定 差异 ,因此 借助 直接 对 比 法 的 定性 
判断 结果 发 现 , MixSIAR 模型 解析 坡度 24" 样 地 柠 条 
绵 鸡 儿 主 要 水 分 来 源 利用 比例 的 适用 性 高 于 Tso- 
Source 模 型 。 曾 祥 明 等 "通过 直接 相关 法 判断 植物 
主要 吸水 源 来 衡量 IsoSource 模 型 和 MixSIAR 模型 
的 适用 性 ,结果 也 表明 利用 MixSIAR 模 型 计算 植物 
对 各 水 源 的 利用 比例 适用 性 高 于 IsoSourece 模型 。 
同时 根据 模型 性 能 评价 结果 ,IsoSource 模型 和 Mix- 
STAR 模型 分 别 计算 坡度 6° .10? 样 地 柠 条 绵 鸡 儿 水 
分 来 源 时 具有 更 好 的 植物 水 分 溯源 效果 ;而 在 计算 
坡度 16° .24? 样 地 柠 条 绵 鸡 儿 水 分 来 源 时 ,IsoSource 
模型 的 准确 性 低 于 MixSIAR 模型 。 因 此 本 研究 认为 
IsoSource 模 型 更 适合 解析 较 小 坡度 (6" 和 16?) 样 地 


柠 条 绢 鸡 儿 的 水 分 利用 策略 ; 而 MixSIAR 模型 解析 
较 大 坡度 (10" 和 24?) 样 地 柠 条 绢 鸡 儿 水 分 利用 策略 
的 适用 性 更 好 。Wang 等 利用 RMSE、G、R、MaxE 
和 MinE 模 型 性 能 评价 指标 对 TsoSource 模 型 和 3 种 
贝 叶 斯 混合 模型 (SIAR 、MixSIR 、MixSIAR ) 的 水 源 分 
配 性 能 进行 了 对 比 评价 ,同样 发 现 IsoSource 模 型 的 
预测 虽然 准确 ,但 该 模型 的 性 能 不 稳定 ,而 SIAR 模 
型 和 MixSIAR 模型 比 MixSIR 模型 具有 更 好 的 植物 
水 分 溯源 效果 ,可 能 是 由 于 后 验 分 布 仿真 算法 的 不 
同 造成 了 这 种 差异 。 

IsoSource 模型 被 广泛 运用 于 解析 植物 水 分 来 
源 ,但 该 模型 只 是 基于 稳定 同位 素质 量 守恒 定律 ， 
未 明确 纳入 水 源 中 同位 素 组 成 的 不 确定 性 ”。 而 
MixSIAR 模型 是 基于 MixSIR 模型 和 SIAR 模型 的 优 
点 所 做 的 改进 ,通过 考虑 源 值 .分 类 和 连续 协 变量 
和 先 验 信息 的 不 确定 性 来 改进 更 简单 的 线性 混合 
模型 ,以 提高 模型 结果 的 准确 性 ""。 同 时 MixSIAR 
模型 提供 了 一 个 最 优 解 , 而 不 是 可 行 解 的 范围 1。 
本 研究 中 坡度 6°、10° 、16° 样 地 柠 条 绵 鸡 儿 木 质 部 水 
及 不 同 土 层 深度 土壤 水 同位 素 组 成 的 空间 差异 性 
较 小 ,即使 IsoSource 模 型 计算 时 未 考虑 植物 木质 部 
水 及 不 同 土 层 深度 土壤 水 同位 素 组 成 的 空间 差异 
性 ,该 模型 计算 结果 的 可 信 度 也 较 高 ,因此 两 种 模 
型 均 适 用 有 旦 计算 结果 的 不 确定 性 较 小 ;而 坡度 24° 
样 地 柠 条 绵 鸡 儿 木 质 部 水 及 不 同 土 层 深度 土壤 水 
同位 素 组 成 的 空间 差异 性 较 大 ,使 得 IsoSource 模 型 
在 计算 坡度 24° 样 地 柠 条 绵 鸡 儿 主 要 水 分 来 源 利用 
比例 时 标准 差 较 大 ,增加 了 该 模型 计算 坡度 24° 样 
地 柠 条 绵 鸡 儿 主 要 水 分 来 源 利 用 比例 的 不 确定 
性 。MixSIAR 模型 解决 了 IsoSource 模 型 计算 时 未 考 
虑 植物 木质 部 水 和 各 潜在 水 源 稳定 同位 素 组 成 空 
间 差 异性 的 问题 ,其 计算 结果 的 标准 差 显 车 低 于 
IsoSource 模 型 计算 结果 的 标准 差 ,有效 降低 了 该 模 
型 计算 结果 的 不 确定 性 , 即 MixSIAR 模型 的 计算 结 
果 能 更 精确 的 给 出 坡度 24° 样 地 柠 条 绵 鸡 儿 主 要 的 
水 分 来 源 以 及 对 不 同 土 层 深度 土壤 水 的 利用 比例 ， 
这 与 张 宇 等 中 的 研究 结果 一 臻 。 李 楠 玉 研究 也 发 
现 ,SIAR 模型 计算 的 各 水 源 贡 献 率 的 标准 差 明 显 低 
于 IsoSource 模型 计算 结果 的 标准 差 。 本 研究 进行 
了 IsoSource 模型 和 MixSIAR 模型 植物 水 分 溯源 效 
果 的 对 比 研 究 ,并 选择 出 了 更 适合 解析 干旱 区 不 同 
坡度 样 地 植物 水 分 来 源 的 溯源 模型 ,解答 了 研究 者 
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在 研究 干旱 区 植物 水 分 来 源 时 ,应 该 选择 哪 种 济源 
模型 的 困惑 。 


5 


结论 


=| 


(1) 柠 条 绵 鸡 儿 对 不 同 土 层 深度 土壤 水 的 利用 


存在 明显 的 季节 性 差异 ,生长 季 初 期 , 随 着 坡度 的 
增 大 , 柠 条 绵 鸡 儿 对 主要 水 分 来 源 的 利用 率 呈 现 出 
先 增 后 减 趋势 ;生长 季 中 期 , 森 条 绢 鸡 儿 主要 水 分 
来 源 随 坡度 的 增 大 由 深 变 浅 ;而 生长 季 后 期 , 柠 条 
绢 鸡 儿 主要 水 分 来 源 随 坡度 增 大 由 浅 变 深 , 呈 相反 
的 变化 趋势 。 


(2) IsoSource 模型 更 适合 解析 较 小 坡度 (6" 和 


16° ) 样 地 柠 条 绢 鸡 儿 的 水 分 利用 策略 ;而 MixSIAR 
模型 解析 较 大 坡度 (10" 和 24?) 样 地 柠 条 绵 鸡 儿 水 分 
利用 策略 的 适用 性 更 好 。 
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A HE: 宁夏 河东 沙 地 不 同 坡度 柠 条 绵 鸡 儿 (Caragana korshinskii) 水 分 利用 策略 差异 


Differences in water use strategies of Caragana korshinskii at different slopes in 
the east sandy land of the Yellow River in Ningxia 


GAO Yang', HAN Lei™, LIU Lili, WANG Nana’, PENG Ling’, 
ZHOU Peng*, ZHAN Xiuli 


(1. School of Geography and Planning, Ningxia University, Yinchuan 750021, Ningxia, China; 2. Joint International Research 
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3. Ningxia Key Laboratory of Resource Assessment and Environment Regulation in Arid Region, Yinchuan 750021, Ningxia, 


China; 4. School of Agriculture, Ningxia University, Yinchuan 750021, Ningxia, China) 


Abstract: We investigated the differences and applicability of IsoSource and MixSIAR models in resolving the 
water use strategies of Caragana korshinskii in the east sandy land of the Yellow River in Ningxia Hui 
Autonomous Region, China. Under the context of global climate change, we used IsoSource and MixSIAR 
models to analyze the use of each potential water source by Caragana korshinskii on different slopes (6°, 10°, 
16°, and 24°) during different periods of the growing season. We then applied hydrogen- oxygen stable isotope 
technique and direct comparison method to evaluate the effectiveness of both models for plant water traceability. 
The results showed seasonal differences in the use of soil water by Caragana korshinskii at different soil depths. 
With the increasing slope, the use of soil water by Caragana korshinskii in the middle layer first increased and 
then decreased at the beginning of the growing season. In the middle of the growing season, the main source of 
water for Caragana korshinskii shifted from deep to shallow soils. Later in the growing season, the main source 
of water for Caragana korshinskii shifted from shallow to deep soils. On the basis of the qualitative results of the 
direct comparison method, both models were applicable for calculating the utilization rate of the main water 
source of Caragana korshinskii in the small sites on 6°, 10°, and 16° samples. The MixSIAR model was more 
reliable than the IsoSource model in calculating the main water source of Caragana korshinskii and its 
contribution rate on the large slope samples. Therefore, the IsoSource model was suitable for analyzing the water 
use strategies of Caragana korshinskii on small slope (6° and 16°) samples, whereas the MixSIAR model was 
suitable for analyzing the water use strategies of Caragana korshinskii on large slope (10° and 24°) samples. The 
results of this study provide a scientific reference for the selection of methods to identify plant water sources in 
arid zones in China. 

Key words: IsoSource model; MixSIAR model; stable isotope; water use strategy; east sandy land of the Yel- 


low River in Ningxia 


